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Abstract. Core-collapsed supernovae (CCSNe) have been considered to be 
one of sources of dust in the universe. What kind and how much mass of dust are 
formed in the ejecta and are injected into the interstellar medium (ISM) depend 
on the type of CCSNe, through the difference in the thickness (mass) of outer 
envelope. In this review, after summarizing the existing results of observations 
on dust formation in CCSNe, we investigate formation of dust in the ejecta and 
its evolution in the supernova remnants (SNRs) of Type II-P and Type lib SNe. 
Then, the time evolution of thermal emission from dust in the SNR of Type lib 
SN is demonstrated and compared with the observation of Cas A. We find that 
the total dust mass formed in the ejecta does not so much depend on the type; 
~ 0.3 - 0.7Mq in Type II-P SNe and - O.lSAfg in Type lib SN. However the 
size of dust sensitively depends on the type, being affected by the difference in 
the gas density in the ejecta: the dust mass is dominated by grains with radii 
larger than 0.03 /zm in Type II-P, and less than 0.006 /im in Type lib, which 
decides the fate of dust in the SNR. The surviving dust mass is ^ 0.04 — 0.2Mq in 
the SNRs of Type II-P SNe for the ambient hydrogen density of nn = 10.0 — 1.0 
cm~'^, while almost all dust grains are destroyed in the SNR of Type lib. The 
spectral energy distribution (SED) of thermal emission from dust in SNR well 
reflects the evolution of dust grains in SNR through erosion by sputtering and 
stochastic heating. The observed SED of Cas A SNR is reasonably reproduced 
by the model of dust formation and evolution for Type lib SN. 



1. Introduction 

Dust is one of important ingredients in space to control the physical and chem- 
ical conditions in interstellar medium (ISM) and the formation process of stars 
in molecular clouds via the interactions with radiation and gas as well as the 
formation of H2 molecule on the surface. In astrophysical environments, dust 
grains form in a cooling gas outflowing from star to interstellar space such as in 
the stellar winds from asymptotic giant branch (AGB) stars and in the ejecta 
of novae and supernovae (SNe) where gas density is high enough to proceed the 
formation and growth of seed nuclei through collisions of relevant gas species. 
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Formation of dust in the ejecta of SNe had been invoked to explain the origin 
of interstellar dust (Cernuschi, Marsicano, & Codina 1967), and the microscopic 
analysis of individual grains in meteorites has identified the dust grains formed in 
SNe from their isotopic compositions (see Clayton & Nittler 2004 for a review). 
Dust formation in the ejecta of SNe had been directly observed for the first 
time in SN 1987A (see McCray 1993 for a review), after then dust formation 
in the ejecta was reported in a handful of core-collapsed SNe (CCSNe). Based 
on the observations of nearby SNe, the mass of dust formed in the ejecta has 
been considered to be less than 10~^Mq per SN, which is more than two to 
three orders of magnitude smaller than the dust mass predicted by theoretical 
investigations (Kozasa, Hasegawa, &: Nomoto 1989, 1991; Todini & Ferrara 2001; 
Nozawa et al. 2003). Furthermore the observed dust mass is too small to explain 
the amount of dust observed in the host galaxies of quasars at redshift of z > 5 
(e.g. Bertoldi et al. 2003) where CCSNe are considered to be a major source of 
dust because the cosmic age is too young for AGE stars to supply the dust. So 
how much mass of dust forms in the ejecta has been still controversial. 

In the ejecta of SNe, dust formation proceeds within the He-core where 
condensible elements are more abundant and the gas density is higher enough 
due to the smaller expansion velocity than in the outer envelope. How much 
mass and what kind of dust are injected from a SN into ISM depend not only on 
the formation in the ejecta but also on the destruction in the supernova remnant 
(SNR) where dust grains are injected into and eroded by sputtering in the hot 
gas swept up by the reverse and forward shocks produced by the interaction 
of ejecta with interstellar/circumstellar medium. Both processes depend on the 
type of SNe through the thickness (mass) of outer envelope: With the same 
kinetic energy of explosion, the thicker outer envelope makes the expansion 
velocity slower and results in the gas density within He-core being high enough 
to form large-sized dust. Furthermore, the lower gas density in the shocked 
region caused by the delayed arrival of the reverse shock at the dust forming 
region decreases the erosion rate of dust by sputtering. 

In this review we show how the formation and destruction of dust in SNe 
depend on the type, comparing the results of calculations for Type II-P and 
Type lib SNe; Type II-P SN is a typical type among CCSNe, and the well 
studied Cas A was recently identified as Type lib SN (Krause et al. 2008). 
Then, we demonstrate the time evolution of thermal emission from dust residing 
in SNR of Type lib SN, based on the results of calculations of dust formation 
and destruction. In § 2, we summarize the observations of dust formation in 
SNe after briefly introducing the classification of SNe. We present the results 
of calculations of dust formation in the ejecta of Type II-P and lib SNe in § 3. 
Focusing on Type lib SN, in § 4 we calculate the evolution of dust in the SNRs 
and the time evolution of thermal emission from dust, and compare with the 
observation of Cas A. Summary and concluding remarks are presented in § 5. 

2. Observations of dust formation in supernovae 

In this section we summarize the observations of dust formation in SNe. First we 
briefly introduce the classification of SNe for your information, and then review 
the observations of dust forming SNe. Massive stars whose main sequence mass 
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is larger than 8 Mq end up their hves as CCSNe, but the details depend on the 
initial metallicity (Heger et al. 2003 for details). Here we shall confine to SNe 
evolved from stars with solar metallicity. 

2.1. Classification of supernovae 

Supernova is the most energetic explosion in the universe. Physically the su- 
pernovae are classified by the explosion mechanism; Core-collapsed supernovae 
which are triggered by the gravitational collapse of the central core, and ther- 
monuclear supernovae in which the explosion energy is supplied by explosive 
thermonuclear burning (e.g. Hillebrandt & Niemeyer 2000). 

Observationally SNe are classified by the spectral features at the early phase 
and the behavior of the light curve (e.g. Wheeler & Harkness 1990; Filippenko 
1997; Leibundgut 2008). According to the presence of the H feature, SNe are 
divided into Type I and Type II; Type I without H, and Type II with H. Type I 
SNe are classified by the strength of the Si feature (Si II absorption at 6150 A); 
Type la SNe showing the strong Si absorption feature are thermonuclear SNe 
whose progenitors are white dwarfs in binary systems. H-deficient Type I SNe 
with weak Si feature are subdivided by the content of He. He-rich is Type lb 
and He-poor is Ic. The progenitors are Wolf-Rayet stars; H-envelope, and both 
H- and He-envelopes are removed before explosion for Type lb SN and Type Ic 
SN, respectively, by stellar wind and/or interaction in binary system. 

Type II SNe, whose progenitors are red supergiants, are divided by the 
content of H; H-poor are lib SNe in which almost all of the H-envelope is re- 
moved during the evolution. H-rich Type II SNe are classified by the behavior 
of the light curves; the light curve linearly decays in Type II-L SNe, and have 
a plateau in Type II-P SNe. Type II SNe showing narrow emission lines origi- 
nating in circumstellar medium are labelled as Iln. Type II, lb and Ic SNe are 
core-collapsed supernovae. 

2.2. Observations of dust formation in core collapsed supernovae 

Dust formation in the ejecta has been observed so far in the ejecta of CCSNe 
except for Type Ic SNe. No evidence of dust formation is reported in the ejecta 
of Type la SNe. Table 1 summarizes the observations of dust formation in the 
ejecta of CCSNe. Formation of dust in the ejecta had been confirmed observa- 
tionally for the first time in SN 1987A which is a peculiar Type II SN because 
the progenitor is not a red supergiant but a blue supergiant (West el al. 1987). 

The evidences of dust formation in the ejecta of SN 1987A (see McCray 
1993; Wooden 1997; the references in Table 1) are the following; (1) the en- 
hancement of infrared fiux at ~ 10 /im starting from day 450 followed by the 
decline of U to M bolometric luminosity from ~ day 500, and (2) the appearance 
of blue-shifted profiles of [Mg I], [O I], [C I] lines starting from ~ day 530. The 
behavior of bolometric luminosity obtained by combining the observed U to M, 
infrared (thermal emission from dust). X-ray and 7-ray luminosities satisfies 
the energy budget expected from the decay of radioactive element ^^Co. The 
appearance of blue-shifted lines was interpreted as the result of the attenua- 
tion of emission from the far (receding) side of the ejecta caused by dust grains 
formed in clumps within the ejecta, and in addition the intensities of [Mg I], [O 
I] and [Si I] lines concurrently decreased (Lucy et al. 1991). It should be noted 
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that, in prelude to the dust formation, formation of CO and SiO molecules was 
observed as early as day 112 and 165, respectively. The emission of SiO was no 
longer recognized at day 517 around which dust formation was confirmed by the 
observation of blue-shifted lines (Roche et al. 1991). The dust mass estimated 
from the observation was at least 1 x 10~^Mq at day 775, and no information 
on the dust composition was available since the thermal emission was well fitted 
by a grey body (Wooden et al. 1993). 

Taking the appearances of the blue-shifted lines, infrared excess and/or 
decline of optical luminosity as the diagnostics, so far dust formation in the 
ejecta has been reported in four SNe excluding SN 1987A (see the references 
in Table 1); three are Type II-P SNe which is the most common type among 
CCSNe (e.g. Smartt et al. 2009), and one is Type lb. The evidence of dust 
formation in SN 2005af is based only on the mid-infrared (MIR) excess observed 
by Spitzer, but the detection of CO and SiO molecules supports the formation 
since CO and/or SiO molecules are detected prior to dust formation in Type 
II-P SNe. The molecules could play an important role in formation of dust 
because CO acts as a coolant in the ejecta (e.g. Liu &: Dalgarno 1995), and SiO 
is considered to be a precursor of dust grains (Kozasa et al. 1989). The onset 
of dust formation ranges from ~ 300 to ~ 600 days after the explosion in Type 
II-P SNe. In SN 19901, which is Type lb, the onset seems to be earlier, as can 
be expected, because the escape of 7-rays with less effective deposition causes 
the gas in the ejecta to cool down faster without an H-envelope. The dust mass 
derived from the observations is less than 1O~^M0, although the observation of 
SN 1999em suggested that dust much more than IO-^'Mq was produced in the 
ejecta (Elmhamdi et al. 2003). 

Recently formation of dust in cool dense shells (CDSs) generated by the 
interaction of ejecta with a circumstellar medium (CSM) has been reported in 
three SNe summarized in Table 2 (see the references in Table 2); two are Type 
Iln and one is Type lb, where the presence of dense circumstellar envelopes is 
confirmed by the X-ray observations (e.g. Smith et al. 2008a). The evidences 
come from the NIR/MIR excess and the concurrent appearance of the blueshift of 
the narrow lines in CSM and/or the intermediate-width components originating 
in shocked gas. Onset of dust formation depends on when the ejecta encounters 
with dense CSM. The early formation of dust at day ~ 50 in SN 2006jc is 
consistent with the LBV-like outburst ~ 2 yr prior to the explosion. The dust 
mass estimates are > 2 x IO'^Mq in SN 1998S, and 3 x 10"'' by day 230 in 
SN 2006jc, assuming carbon dust. It should be noted that in addition to the 
formation of dust in CDSs, we cannot deny the possibility that dust also forms in 
the expanding ejecta; see Smith et al. (2008b) for SN 2005ip; Sakon et al. (2009), 
Tominaga et al. (2008), and Nozawa et al. (2008) for SN 2006jc. Although the 
frequency of Type Iln SNe is rare among CCSNe (Gal- Yam et al. 2007; Smartt 
et al. 2009), these observations on dust formation in CDS provide a new window 
to investigate the formation process of dust in astrophysical environments. 

Nowadays, based on the observations, the mass of dust formed in the ejecta 
has been claimed to be less than 10~^Mq and to be too small to contribute 
to the inventory of dust in our Galaxy. However, we should keep in mind that 
the conclusion is not necessarily definite, based on the limited number of the 
observations and the assumption that thermal radiation from dust is optically 



Table 1. Formation of dust in the ejecta. References; 1) Lucy et al. (1991), 2) Meikle et al. (1993), 3) Roche et al. (1989), 
4) Whitelock et al. (1989), 5) Suntzeff & Bouchct (1990), 6) Lucy et al. (1989), 7) Catchpole & Glass (1987), 8) Aitken et al. 
(1988), 9) Ercolano et al. (2007), 10) Elmhamdi et al. (2004), 11) Elmhamdi et al. (2003), 12) Spyromilio et al. (2001), 13) 
Sugerman et al. (2006), 14) Hendry et al. (2005), 15) Meikle et al. (2007), 16) Kotak (2008), 17) Kotak ct al. (2006). Note; 
The optical depth of the dust core derived from [O I] line profile Td ^ 10 at day 510 in SN 1999em (Elmhamdi. ct al. 2003), 
which is much larger than Td ~ 1 in SN 1987A (Lucy et al. 1991). 



name 


Type 


onset of dust 


IR excess 


decline of 


blue-shifted 


molecules 








formation tc[day] 


(NIR/MIR) 


opt. lum. 


lines 




[Mo] 


SN 1987A 


Il-pec. 


~ 350-530 i)'2) 


MIR 3) 


yes 4)'5) 


[Mg I],[0 I] etc. 6) 


CO ^), SiO 8) 


7.5 X 10-4 9) 


SN 19901 


lb 


- 250 




yes 10) 


[0 I], [Call] 10) 






SN 1999em 


II-P 


465< tc <510 11) 




yes 11) 


[0 I],[Mg I], etc. 11) 


CO 12) 


>10-4 11) 


SN 2003gd 


II P 


250< tc <493 13) 


MIR 13) 


yes 13) 


[0 I], H„ 13)'14) 


CO,SiO ? 15) 


4 X 10-5 15) 


SN 2005af 


II P 


214< tc <571 1^) 


MIR 16) 






CO,SiO 1^) 


-4x10-4 16) 



Table 2. Formation of dust in a cool dense shell produced by the interaction with dense circumstellar medium. References; 1) 
Pozzo et al. (2004), 2) Gerardy et al. (2000), 3) Smith et al. (2008b), 4) Fox et al. (2009), 5) Smith et al. (2008a), 6)Mattila 
et al. (2008). Note: For SN 2005ip, Fox et al. (2009) suggested the early formation of dust at day ~ 50, but the spectroscopic 
observation (Smith ct al. 2008b) suggested the later formation. It should be noted that Smith ct al. (2008b) have suggested 
formation of dust in the ejecta at earlier time. 



name 


Type 


onset of dust 


IR excess 


decline of 


blue-shifted 


molecules 


Md 




formation tc [day] 


(NIR/MIR) 


optical lum. 


lines 




[Me] 


SN 1998S 


Iln 


- 250 1) 


NIR 1) 




H„, He I 1) 


CO 2) 


> 2 x 10-3 1) 


SN 2005ip 


Iln 


- 730 3) 


NIR 3^4) 




He I 3) 






SN 2006jc 


lb 


~ 50 5'6) 


NIR-MIR '^.e) 


yes 


He I 5'6) 




3 X 10-4 6) 
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thin in some cases; as pointed out by Meikle et al. (2007), the dust may reside in 
optically thick clumps in the ejecta. Also, the NIR-MIR observations could miss 
the cool dust, since dust grains might cool down quickly after the formation as 
is demonstrated by Nozawa et al. (2008). Therefore, the monitorings of various 
types of SNe after the explosions, covering the temporal and wavelength ranges 
as wide as possible, are necessary to reveal the mass of dust formed in the 
ejecta and its dependence on the type of SNe. In addition the sophisticated 
radiative transfer calculations such as those carried out by Sugerman et al. 
(2006) and Ercolano et al. (2007) are inevitable to derive the dust mass from the 
observations relevant to dust formation, taking into account the more realistic 
spatial distribution and chemical composition of dust in the ejecta. 

3. Formation of dust in the ejecta 

In this section we investigate how the formation process of dust in the ejecta 
depends on the type of SNe, focusing on Type II-P and Type lib SNe whose 
progenitors are evolved from massive stars with solar metallicity. The calculation 
of dust formation is based on the theory of nucleation and growth developed by 
Nozawa et al (2003), taking into account chemical reaction at the condensation 
by considering that the key species, defined as a gas species with the least 
collisional frequency among the reactants, controls the kinetics of nucleation and 
grain growth (Kozasa & Hasegawa 1987, Hasegawa & Kozasa 1988) Q. Given 
the time evolution of gas density and temperature together with the elemental 
composition, we can determine when, where, what kind, what size, and how 
much mass of dust condenses in the ejecta. 

3.1. Models of supernovae and dust formation calculation 

Formation process of dust grains in the ejecta of SNe is controlled by the time 
evolution of gas density and temperature as well as the elemental composition 
in the ejecta as demonstrated by Kozasa et al. (1989). Thus, we apply the 
hydrodynamic model of a SN explosion for the time evolution of the gas density 
and the nucleosynthesis calculation for the elemental composition in the ejecta. 
The time evolution of the gas temperature is calculated by solving the radiative 
transfer equation taking into account the energy deposition from radioactive 
elements. Table 3 summarizes the models of SNe used for the calculations. 

In the calculations we consider that dust forms within the He-core because 
the large expansion velocity in the outer envelope causes the density of conden- 
sible elements to be too low to form dust grains. For the elemental composition, 
we do not consider any mixing since the mixing driven by the Rayleigh-Taylor 
instability in the ejecta at the explosion is not microscopic but macroscopic as 
have been revealed by the observation of Gas A (Douvion et al. 1999; Funis et 
al. 2006). We assume that formation of CO and SiO molecules is complete; all 
carbon (silicon) atoms are locked into GO (SiO) molecules in the locations where 
C/0 (Si/0) is less than 1. This implies that G-bearing dust condenses only in 



^ There was a typo in the nucleation rate J in Kozasa & Hasegawa (1987) and Nozawa et al. 
(2003); replace ci with Cj, where ci is the number density of key species. 
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Table 3. Models of SNe; Mpr is the progenitor mass, ii^oxp the kinetic energy 
of the explosion, M(^^Ni) the mass of ^^Ni in the ejecta, Afcjc the mass of the 
ejecta, Mn-env the mass of the hydrogen envelope, Mhc-coic the mass of the 
He-core, and Mmatai the mass of metal inside the He-core. Model A is taken 
from Nomoto et al. (2006), and model B from Umcda and Nomoto (2002). 



model 


Type 








Mojc 




-^-^Hc— core 


A/mctal 








[erg] 




[Mq] 


[Mq] 


[Me] 


[Mq] 


A 


II-P 


15.0 


10^1 


0.134 


12.2 


10.4 


1.80 


0.482 


B 


H-P 


20.0 


10^1 


0.07 


15.7 


10.9 


4.77 


3.35 


C 


Hb 


18.0 


10^1 


0.07 


2.94 


0.08 


2.86 


1.297 



the region of C/0> 1; see Todini & Ferrara (2001) and Bianchi &; Schneider 
(2007) for the calculation of dust formation assuming the uniform mixing of 
elements and considering formation processes of CO and SiO in the ejecta. 

Figure 1 shows the elemental composition within the He-core for a Type 
II-P SN (model A; upper panel) and for a Type lib SN (model C; lower panel). 
From the view point of dust formation, the He-core is roughly divided into four 
layers; the outer C-rich layer, the 0-rich layer, the Si-S-Fe layer and the Fe- 
S-Si layer. What kind of dust grain really condenses depends on the details 
of elemental composition. Thus, in the calculation we simultaneously solve the 
equations of nucleation and grain growth for 19 possible condensates (see Nozawa 
et al. 2003 for details). 

Figure 2 shows the time evolution of the gas density and temperature in 
the ejecta; Left panel for Type II-P SN (model A) and right panel for Type lib 
SN (model C). It should be noted that the gas density in the ejecta of Type 
lib SN is more than three orders of magnitude smaller than that of type II-P 
SN, because the expansion velocity in the ejecta is much higher in Type lib SN 
without thick H-envelope. The time evolution of the gas temperature depends 
on the model through the mass of H-envelope and M(^^Ni) in the ejecta. 

3.2. Results of dust formation calculations 

Figure 3 shows the condensation time of each grain species formed in the ejecta; 
Left for Type II-P SN (model A) and right for Type lib SN (model C). As the 
gas in the ejecta cools down with time, C grains start to condense first in the 
C-rich layer around day 300 after the explosion. Afterwards, in the oxygen-rich 
layer, AI2O3, Mg2Si04, MgSiOa and Si02 condense in this order from day 350 
to day 450. In the Si-S-Fe rich layer, FeS and Si condense around day 450 in 
Type II-P SN and around day 630 in Type lib SN. After then, in the innermost 
Fe-S-Si layer, Si, FeS and Fe condense in Type II-P SN, but Si and FeS in Type 
lib SN. The condensation sequence of MgO and Fe304 in the 0-rich layer is 
different between Type Il-P SN and Type lib SN, depending on the detailed 
density structure and elemental composition in the ejecta. Also small amount 
of FeS and Si condenses from 450 to 500 days in the C-rich layer of Type II-P 
SN, and so does FeS around day 580 in the 0-rich layer of Type lib SN. 
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Figure 1. Elemental composition within He-core: Upper panel for Type 
II-P SN (model A) and lower panel for Type lib SN (model C). Note that 
the position in the ejecta is indicated by the mass coordinate that is defined 
as the enclosed mass from the center. 
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Figure 3. Condensation times of dust species formed in the cjecta; Left 
for Type II-P SN (model A) and right for Type lib SN (model C). The 
condensation time is defined as the time at which the nucleation rate reaches 
the maximum. The condensation times of dust species in the C and 0-rich 
layers are a little earlier in Type II-P than in Type lib, except for FeS and 
Si in C-rich layer of Type Il-P SN, and FeS in Type lib 



Figure 4 shows the average radius of each dust species formed in the ejecta. 
In Type lib SN (right panel), the average radii of all grain species are smaller 
than 0.01 /xm due to the low gas density in the ejecta. On the other hand, the 
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1.5 2 2.5 3 2 3 4 

mass coordinate; (^q) mass coordinate; (Mg) 

Figure 4. Average radii of dust species formed in the ejecta: Left for Type 
II-P SN (model A) and right for Type lib SN (model C) 




grain radius {/^m) grain radius (/.im) 



Figure 5. Cumulative size distributions of dust species by mass summed up 
over the ejecta. The mass is dominated by dust grains with radii larger than 
0.03 /im for Type II-P SN (model A; left), but with radii less than 0.006 ^im 
for Type lib SN (model C; right) 



average radii of dust grains are larger than 0.01 fim in Type II-P (left panel), 
except for Si and FeS in the C-rich layer, and AI2O3, MgO and Fe304 in the 
0-rich layer. The difference in the size of dust grains between Type Il-P SN and 
Type lib SN is much more clearly represented by the cumulative size distribution 
of dust by mass summed up over the ejecta, which is given in Figure 5. Except 
for AI2O3, the dust mass is dominated by grains with radii larger than 0.03/im 
in Type II-P SN (model A; left panel), while the dust mass of all species is 
dominated by the grains whose radii are less than 0.006 //m in Type lib (model 
C; right panel). This leads to the big difference in the fate of dust grains in 
SNRs as is presented in the next section. 

Table 4 summarizes the mass of dust formed in the ejecta of Type II-P 
SNe (model A and B) and Type lib SN (model C). The masses of C, Mg- 
silicates (Mg2Si04 and MgSiOs), Si02, Si and FeS do not depend so much on 
the model, while masses of AI2O3, MgO, Fe304 and Fe are different from model 
to model, depending on the detailed elemental composition and density structure 
in the ejecta. Generally the total dust mass increases with increasing progenitor 
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Table 4. Dust mass formed in the the ejecta 



model 


A 




B 




C 




Type 


n-p 




II-P 




lib 




MpAMq] 


15 




20 




18 




dust species 














C 


4.30 X 10" 


2 


4.00 X 10- 


2 


7.08 X 10- 


2 


AI2O3 


1.76 X 10" 


3 


4.33 X 10- 


2 


6.19 X 10- 


5 


Mg2Si04 


2.95 X 10" 


2 


0.133 




1.74 X 10- 


2 


MgSiOa 


1.91 X 10" 


3 


5.68 X 10- 


3 


5.46 X 10- 


2 


MgO 


3.19 X 10" 


7 


0.159 




2.36 X 10- 


3 


Si02 


2.06 X 10"^ 


2 


8.21 X 10- 


2 


1.57 X 10- 


2 


Fe3 04 


5.60 X 10- 


5 


1.55 X 10- 


3 






FeS 


3.79 X 10- 


2 


5.66 X 10- 


2 


1.47 X 10- 


3 


Si 


6.65 X 10- 


2 


9.04 X 10- 


2 


5.07 X 10- 


3 


Fe 


0.134 




6.46 X 10- 


2 






total mass 


0.327 




0.676 




0.167 




Md/A/mctal 


0.695 




0.202 




0.129 






1.5 3 2.5 3 1.5 3 3.5 3 3.5 4 

mass coordinate; M^. (Mg) mass coordinate; {^q} 



Figure 6. Optical depth of the ejecta at 1000 days after explosion as a 
function of mass coordinate; Left for Type 11-P SN (model A) and right for 
Type lib SN (model C). 

mass. The total dust mass ranging from ~ 0.17 to 0.68 Mq is two to three 
orders of magnitude larger than the mass inferred from the observations. The 
condensation efficiencies defined by Md/Mmetai ranges from 0.13 (model C) to 
0.7 (model A). The condensation efficiency of model A is significantly higher 
than in the other models, because the metal mass within the He-core is smaller 
due to the thin 0-rich layer (see Table 1 and Figure 1). Thus, being different 
from the total dust mass, the condensation efficiency strongly depends on the 
detailed density structure and elemental composition in the ejecta. 

Finally we present the optical depth of the ejecta at day 1000 after the 
explosion in Figure 6. The optical depth of the ejecta after dust formation 
decreases proportional to since the ejecta expands homologously. Thus, 
around day 500 after the onset of dust formation, the ejecta is optically thick 
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even at 24 /xm except for the outer ^ 0.3 Mq (1.5 Mq) region for model A 
(model C). This result suggests that it could be inappropriate to estimate dust 
mass from the thermal emission by assuming the ejecta to be optically thin at 
MIR. It takes about 800 (60) yr at 0.5 /im for the entire ejecta to be optically 
thin, and about 50 (10) yr at 24 ^um for model A (model C). 

4. Evolution of dust in supernova remnant 

Dust grains formed in the ejecta suffer from erosion and destruction by sput- 
tering in the hot gas between the reverse and forward shocks produced by the 
interaction of the ejecta with the surrounding medium. This process determines 
how much mass and what kind of dust can survive and be injected into ISM. In 
this section we demonstrate the evolution of dust in SNRs, focusing on type lib 
SN (model C). Then, we calculate the time evolution of thermal emission from 
dust in the SNR and compare with the observation of Cas A. 

4.1. Motion and destruction of dust in SNRs 

At the time when the ejecta encounters the ISM surrounding the SN, the for- 
ward shock propagating into the ISM and the reverse shock into the ejecta are 
generated (e.g. Truelove and McKee 1999), and the gas swept up by the shocks 
is heated up and ionized. The dust grains hit by the reverse shock are injected 
into the shocked gas with a velocity relative to the gas; gas is decelerated by 
the shock, while dust grains decoupling from gas penetrate into the hot gas with 
the same velocity just before the encounter. Then dust grains are decelerated 
and eroded by collision with gas. The deceleration of dust by the gas drag is in- 
versely proportional to pdOj and the erosion rate of dust by sputtering is almost 
independent of a (e.g. Dwek &: Arendt 1992), where a and pd is the radius and 
the bulk density of dust, respectively. As a result, the small sized/low density 
dust is easily trapped into gas, and is eroded by non-kinetic (thermal) sputter- 
ing, while larger sized/high density dust moves through the gas, suffering kinetic 
sputtering. 

We calculate the motion and destruction of dust in the SNR, adapting the 
model for the SN explosion as the initial conditions for gas velocity and density, 
and applying the results of dust formation calculations as the initial spatial and 
size distributions of dust grains in the ejecta. The method of calculation and 
the underlying assumptions are the same as those by Nozawa et al. (2007). In 
the calculations, we assume that the ejecta expands into a uniform medium with 
solar metallicity. We consider the two cases for the ambient hydrogen number 
density; nn =1.0 and 10.0 cm We assume the gas temperature in the ambient 
medium Tg amb = 10^ K, and that the ejecta hits the ISM at tone = 10 yr after 
the explosion. Note that the assumptions on Tg amb and tcnc do not significantly 
affect the result of calculations. The sputtering yield of each grain species is 
taken from Nozawa et al. (2006). 

The left panel of Figure 7 shows the time evolution of gas density and 
temperature in the SNR of type lib SN (Type lib SNR) for nn = 10 cm~^; the 
positions of the reverse (forward) shock at given times are indicated by upward 
(downward) arrows. The temperature of gas hit by the shocks quickly rises up 
to 10^ to 10^ K, and is kept to be > 10^ K until the cool dense shell forms behind 
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the forward shock att ^ 10^ yr after the explosion. The reverse shock encounters 
the C-rich layer at ~ 80 yr, the 0-rich layer at ~ 440 yr, the Si-S~Fe rich layer 
at ^ 800 yr, and the Fe-Si-S rich layer at ^ 880 yr after the explosion, and the 
ejecta are completely swept up by the reverse shock in ~ 1400 yr. The right- 
upper panel shows the trajectories of dust species with radii of 0.001 /.tm (dashed 
curves) and 0.003 fim (solid curves) together with the positions of the forward 
and reverse shocks and the outer boundary of He-core. The time evolution of 
the radius of each grain species is depicted in the right-lower panel. As can be 
seen, grains with radii of 0.003 //m are quickly trapped into the hot gas behind 
the reverse shocks and are completely destroyed by thermal sputtering within 
~ 3 X 10^ yr. How long dust grains remain in SNR depends not only on the 
initial position but also on the chemical composition. 

The time evolution of the mass of each dust species in the hot gas of Type 
lib SNR is displayed in Figure 8; Left for nn = 1-0 cm~^ , and right for nu = 10.0 
cm~^. Except for C grains for nn = 1.0 cm~^, all dust grains are destroyed in 
t < 3x 10^ (4x 10^) yr for rin = 1.0 (10.0) cm~^. Although the time necessary for 
the reverse shock to reach the dust forming layer is almost two times longer for 
nn = 1.0 cm~'^ than for njj = 10.0 cm"'^, it takes much longer for dust grains 
to be destroyed for nn = 1.0 cm~^ because the lower density in the hot gas 
caused by the delayed encounter reduces the erosion rate of dust by sputtering 
significantly. It should be noted that C grains dominate the dust mass almost 
over all times for nn = 1.0 cm~^, while for nn = 10.0 cm~^ the mass of C grains 
quickly decreases at t > 300 yr, and then MgSiOs and Si02 grains dominate the 
dust mass in the SNR. Thus, the ambient gas density strongly affects not only 
the mass but also the composition of dust grains remaining in the SNR. 

Table 5 presents the surviving dust mass Msurv at t = 10^ yr for the models. 
Although Msurv decreases with increasing nn, how much mass of dust survives 
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Figure 8. Time evolution of dust mass in the shocked gas of the Type lib 
SNR: Left for nn = 1-0 cm^'^ and right for nn = 10.0 cm~'^. 



Table 5. Mass of dust surviving in the SNRs Msurv at t=10 ^ yr and the 
efhciency of survival defined by 77 = Msui-v/Md 



nH[cm ^] 




1.0 




10.0 


model (type) 




] V 


M,u,v[Af0 


] V 


A (Type II-P) 


0.190 


0.581 


0.104 


0.318 


B (Type II-P) 


5.66 X 10- 


2 8.37 X 10-2 


3.91 X 10" 


2 5.78 X 10-2 


C (Type lib ) 


7.83 X 10- 


5 6.07 xlO-'' 


0.00 


0.00 



strongly depends on the type of SNe through the thickness of outer H-envelope 
that influences the size of dust formed in the ejecta as weh as the time when 
the reverse shock reaches the dust forming layers. In the Type lib SNR, dust is 
destroyed quickly partly because the small sized dust grains with radii < 0.006 
fim populate the remnant region. In addition, the earlier arrival of reverse shock 
at the dust-forming layer causes the gas density behind the reverse shock to be 
higher, and enhances the erosion rate of dust grains; for example, the 35 times 
earlier arrival to the C-rich layer compared to the Type II-P SNR (model A) 
results in the 10 times higher gas density behind the reverse shock in the Type 
lib SNR. However, it should be noted that, contrary to the expectation, dust 
grains are more destroyed in model B than in model A, despite that the mass 
of H-envelope is almost same. One of the reasons is that MgO grains with a 
sputtering rate somewhat larger than others grains (Nozawa et al. 2006) are 
more abundant in model B (see Table 4) , and another is that the sizes of Fe and 
Si grains in model A are significantly larger than those in model B. Thus, the 
survival of dust grains in SNR depends on not only the thickness of H-envelopes 
but also on the detailed chemical composition, size distribution, and the amount 
of dust grains formed in the ejecta. 

The efficiency of survival i] defined as Msuw/Mfi is larger than ~ 0.06 in 
the Type II-P SNR, but less than 6 x lO"'' in the Type lib SNR. Dust grains 
surviving at t = 10^ yr are almost free from the destruction by sputtering, and 
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are trapped into the dense cool shell formed behind the forward shock or injected 
into the ISM. Thus, the injection efficiency strongly depends on the type, which 
should be incorporated into chemical evolution models of dust in the ISM. In 
addition, it should be mentioned here that the survival in SNRs and the injection 
into the ISM of newly formed dust depend on the density structure in the ISM 
as well as the inhomogeneity in the ejecta. For the ambient medium whose gas 
density decreases as r~" with distance r, the time evolution of the density of gas 
swept up by the reverse shock is different from that for a constant density as 
shown by Chevalier & Oishi (2003), which affects the evolution of dust in SNRs. 
If dust grains reside in dense clumps (knots) as suggested by the observations 
(e.g. Arendt, Dwek, & Moseley 1999), the dust grains could survive, since the 
reverse shock encountering dense clumps could indeed be radiative so that the 
gas temperature could stay cool enough for dust grains to be rescued from the 
destruction by sputtering. These aspects should be explored to reveal how much 
mass of dust grain can survive in SNRs and be injected into the ISM. 

4.2. Thermal emission from dust in Type lib SNR 

Dust grains formed in the ejecta and injected into the hot plasma between the 
reverse and forward shocks are heated up by collision with gas, and emit thermal 
radiation. The equilibrium temperature of dust grain is determined by balancing 
the coUisional heating with the radiative cooling. However, in a rarefied hot 
plasma such as in SNRs, small sized dust grains undergo stochastic heating which 
affects the emissivity and the resulting spectral energy distribution (SED); see 
Dwek (1986) and Dwek et al. (2008) for details. 




Temperature (K) Wavelength A((tmi) 



Figure 9. Temperature distribution functions of small astronomical silicate 
in the hot plasma with = lOcm"'^ and To = 10^ K (left) and the emissivities 
(right) for given radii, where the radius is in units of fim. 

For example. Figure 9 shows the temperature distribution function and the 
emissivity of astronomical silicate (Draine and Lee 1984) embedded in a hot 
plasma with the electron temperature of Tc = 10^ K and number density of 
Tic = 10.0 cm~'^, which is calculated by a Monte-Carlo method. The fraction 
of incident energy that is deposited in a grain is calculated by fitting the ex- 
perimental data of the stopping ranges of electrons tabulated by Iskef et al. 
(1983), according to the method by Dwek (1987b). With decreasing radius, the 
temperature distribution function becomes broader and the resulting emissivity 
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deviates strongly from the emissivity calculated by the equilibrium temperature; 
the deviation becomes significant for the radius of a < 0.02;um. 

Taking into account the stochastic heating, we demonstrate the time evo- 
lution of thermal emission from dust grains embedded in the shocked gas in 
the Type lib SNR, based on the results of dust formation and evolution cal- 
culations for Type lib SN presented in § 3 and § 4.1. In the calculations, the 
heat capacities of grain species are taken from Takeuchi et al. (2005), and we 
use the same optical constants as those of Hirashita et al. (2005), except for 
amorphous AI2O3 (Begemann et al. 1997 for the wavelength A > 8 /im), Si02 
(Phillip 1985) and Si (Pillar 1985). We consider only collision with electrons 
for heating, because the radii of dust grains formed in the ejecta are so small 
that dust grains injected into the hot gas are quickly trapped and comove with 
the gas, and the heating is dominated by collision with electrons. In order to 
compare with the observations of Cas A, we put the SNR at the distance of 3.4 
kpc. Figure 10 and 11 show the time evolution of SED ranging from 200 to 
1800 yr after the explosion for nH=1.0 and 10.0 cm~^, respectively, where the 
solid (dotted) curve denotes the SED with (without) stochastic heating. The 
solid circles are the flux densities subtracting the synchrotron radiation from the 
observed ones tabulated by Hines et al. (2004). 

The SEDs taking into account the stochastic heating is completely different 
from those with the equilibrium temperatures, especially in the wavelengths 
shorter than 20 ^m. After the revere shock encounters the C-rich layer at ~ 
180 (80) yr for nH=1.0 (10.0) cm~^, the thermal emission from the stochastically 
heated C grains contribute to the SED until the reverse shock encounters the 
0-rich layer. In the case of nn = 1-0 cm~^, the thermal emission from C 
grains dominates the SED over the entire simulation period considered in the 
calculation, because the mass of dust in the SNR is almost dominated by C 
grains as can be seen from Figure 8. The emission features of MgO appearing at 
A ~ 16 i^tm from t = 800 to 1200 yr, Mg-silicates at A 10 nm after t = 1200 yr 
and Si02 at A ~ 20 /xm after t = 1600 yr are not prominent, being overwhelmed 
by thermal emission from C grains. On the other hand, in the case of nn = 10.0 
cm~^, with C grains being more efficiently destroyed by sputtering, the emission 
features from other dust grains get prominent with time, after the reverse shock 
encounters the 0-rich layer at t ~ 400 yr; although the feature of MgO is 
weak, the emission feature of Mg-silicates around 10 ^um becomes noticeable 
at t = 600 yr, after then the 20 fj,m feature of Si02 is prominent at t = 800 
yr. With increasing time, the SED is dominated by the emission features from 
stochastically heated dust grains. 

The result of calculations clearly demonstrates that the time evolution of the 
SED is very sensitive to the evolution of dust in SNRs through the destruction 
by sputtering and the stochastic heating, which strongly depends on the density 
in the ambient medium. It should be noted that the calculated SED cannot be 
reproduced by a single dust component with multiple temperatures. 

4.3. Comparison with the observation of Cas A 

Cas A is the nearby Galactic SNR well-studied at various wavelengths and was 
recently identified as Type lib SN from the spectrum of the scattered light 
echo (Krause et al. 2008). The presence of dust formed in the ejecta has been 
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Figure 10. The time evolution of SEDs of Type lib SNR for tih = 1-0 cni-^ 
with (sohd curve) and without (dotted curves) stochastic heating. Closed 
circles are the flux densities subtracting the synchrotron radiation from the 
observed ones of Cas A (Hines et al. 2004) 






Figure 11. The same as Figure 10, but for nn = 10.0 cm 
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Figure 12. Comparison with the observation of Ca A. The symbols denote 
the flux densities of Cas A. The sohd curve is the total flux density, and the 
dotted curves denote the contributions of each grain species. 



confirmed by the infrared and submillimeter observations (Lagage et al. 1996; 
Arendt et al. 1999; Dunne et al. 2003). However, the dust mass derived from 
the observations has been controversial, ranging from < 1 X 10^3 to ~ I.OMq 
(Dwek et al. 1987a; Douvion, Lagage, Pantin 2001; Dunne et al. 2003; Hines 
et al. 2004; Rho et al. 2008). The comparison of the calculated SED with the 
observations of Cas A could provide some insights on dust and its evolution in 
the Cas A SNR, although the models of Type lib SN and the evolution of the 
SNR are not necessarily applicable to the Cas A SNR. 

As shown in Figures 10 and 11, the calculated SEDs for n = 1.0 cm~'^, 
being dominated by thermal emission from small C grains, cannot reproduce the 
observations at the wavelengths A < 20 //m at t < 1000 yr, and are significantly 
smaller than the observed flux densities at A > 20 /xm. On the other hand, for 
nil = 10.0 cm""^, the calculated SED at t = 800 around which the reverse shock 
encounters the Si-rich layer seems to reproduce the observed SED. 

Figure 12 shows the detailed comparison of the SED calculated at t = 800 
yr with the observation and the contribution of each dust grains. Although the 
discrepancy around A ~ 12 ^tm is remarkable and the 20 //m feature of Si02 is 
too sharp, the calculated SED can reasonably reproduce the overall shape of the 
observed SED without any tuning. The SED is dominated by Mg-silicates and 
Si02 at A > 15 /xm, and by C and Mg-silicates at A < 10 //m. The dominant 
dust species (the mass in unit of solar mass) are MgSiOs (1.72 x 10~^), Si02 
(1.08 X 10-2), MgSi04 (6.46 x lO'^), Si (3.62 x 10-^), C (1.6 x 10^^)^ ^nd FeS 
(1.04 X 10^'^). The total dust mass is ~ 0.04 Mq, which is consistent with 
the dust mass of 0.02-0.054 Mq evaluated by Rho et al. (2008), apart from 
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the details. The presence of Si grains in the model indicates that the reverse 
shock already reaches a part of the inner Si-S-Fe layer, which is also not in 
contradiction with the observation of Cas A by Spitzer showing that the deeper 
Si-rich layer has been hit by the reverse shock in the region associated with the 
jet arising from the asymmetric explosion (Ennis et al. 2006). 



5. Summary and concluding remarks 

We calculate dust formation in the ejecta and evolution in the SNRs of Type 
II-P and lib SNe in order to clarify how these precesses depend on the type of 
CCSNe through the difference in the thickness (mass) of the outer envelope. We 
show that the mass of newly formed dust ranging from ~ 0.1 to 0.7 Mq does not 
so much depend on, but the size of dust grains is sensitive to the thickness of 
H-envelope. The radii of dust grains are less than 0.006 nm in Type lib SN with 
Mn-env = 0.08 Mq, wliilc the dust mass is dominated by the grains with radii 
larger than 0.03 /xm in Type II-P SNe with Mn-env ~ WMq. The difference in 
the size of the dust formed in the ejecta plays a crucial role in the evolution of 
dust in SNRs because the smaller sized grains are quickly trapped into the hot 
gas behind the reverse shock and are destroyed by sputtering. The surviving 
dust mass ranges from 0.19 to 3.9 x lO'^M© in the SNR of Type II-P SN for 
the ambient H number density of nn = 1.0-10.0 cm~^, but the dust grains are 
almost completely destroyed in the SNR of Type lib SN. 

The SED of the thermal emission from dust embedded in SNRs can be 
applied as a diagnosis of the evolution of dust, reflecting the destruction by 
sputtering and the stochastic heating. The comparison of the calculated SED of 
Type lib SNRs with the observation of Cas A suggests that the mass of dust in 
the SNR is 0.04 Mq, which implies that the mass of dust formed in the ejecta is 
not less than 1.0 x 10""^ Mq, being different from the dust mass estimated from 
the dust forming CCSNe. In addition, Mg-silicates and Si02 with mass reaching 
up to 3.0 X 10~^Mq dominate the SED at A > 15yum, in contrast to the observed 
dust forming SNe that so far show no signature of silicates. This may indicate 
that these dust grains formed in the 0-rich layer cool down quickly and evade 
being detected by the observations as is discussed by Nozawa et al. (2008). 

It has been claimed that the mass of dust formed per CCSN is less than 
1.0 X 10"'^ Mq and is too small to contribute to the inventory of dust in ISM, 
based on the observations. However, it should bear in mind that the conclusion 
comes from a small number of observations. In order to clarify when how much 
mass and what kind of dust condenses in the ejecta of various types of SNe, it 
is inevitable to monitor SNe with temporal and wavelength coverages as wide 
as possible. In addition, the dedicated radiation transfer calculations such as 
those by Sugerman et al. (2006) and Ercolano et al (2007) are promising in 
evaluating the dust mass compared with the observations, taking into account 
the distribution and chemical composition of dust grains in the ejecta. 
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